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A consistent picture of the Au+Au and D+Au, y/s = 200 A GeV measurements at RHIC ob- 
tained with the PHENIX, STAR, PHOBOS and BRAHMS detectors including both the rapidity 
and transverse momentum spectra was previously developed with the simulation LUCIFER. The 
approach was modeled on the early production of a fluid of pre-hadrons after the completion of 

£ — ' an initial, phase of high energy interactions. The formation of pre-hadrons is discussed here, in a 

perturbative QCD approach as advocated by Kopeliovich, Nemchik and Schmidt. In the second 
phase of LUCIFER, a considerably lower energy hadron-like cascade ensues. Since the dominant 

CNj ' collisions occurring in this latter phase are meson-meson in character while the initial collisions are 

r*. , between baryons, i. e. both involve hadron sized interaction cross-sections, there is good reason to 

suspect that the observed elliptical flow will be produced naturally, and this is indeed found to be 

^yy ' the case. 

PACS numbers: 

I. INTRODUCTION 

RHIC experiments 0, i, i] have identified a transverse momentum asymmetry, designated as 'elliptical flow,' 
and characterized by the variable «2- This quantity has deservedly attracted considerable interest. Early heavy 
i— i , ion experiments at the LBL TEVALAC [J, Q defined flow in terms of the asymmetry of the transverse momentum 
distributions observed for a variety of particles, both baryons and mesons, produced in relativistic heavy ion collisions 
at these relatively low energies. Investigations of flow were subsequently extended to the AGS heavy ion program 
and eventually flow was studied in experiments at the CERN SPS. LBL and even the AGS energies per nucleon were 
of course well below those achievable at RHIC. At lower energies, however, it was possible to examine an event by 
. event phenomenon, and to extract what was generally referred to as directed or sidewards flow Q, separately for 
baryons and mesons. This was of critical importance in understanding the dynamics of relativistic ion collisions at 
, those energies. Much theoretical effort was expended in understanding this phenomenon in terms of hydrodynamics [§} 
employing particular equations of state for high density nuclear matter and setting some set of initial conditions for 
, the hydrodynamics. 

Considerable success, however, was also obtained Q and a description of the sidewards flow using only a cascade 
oriented simulation, also proved to be possible, with its own inherent treatment of the nuclear medium. Here we 
discuss at s 1 / 2 = 200 GeV only the prominent, designated elliptical, flow extracted from a combined history of many 
events. We are also limited to overall particle distributions, necessarily dominated by final stable mesons. 

Flow may be formally introduced by expressing the observed transverse momentum distributions as a function of 
a tranverse azimuthal angle 4>, relative to the angle ipji fixing the inclination of initial reaction plane. 
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Of particular interest at RHIC is the second term in this expansion, i. e. V2, which if non- vanishing indicates an 
elliptical asymmetry. The first term, if non-zero, would signal the presence of directed or sidewards flow, and has been 
estimated in the RHIC detectors [H, [2, 3], the corresponding parameter vt, is found to be small in these experiments. 
Data has been accumulated by three of the RHIC detectors P, 0, [|[ producing experimental studies of the ellipticity 
for different degrees of centrality, and also as a function of p±. A comparison is made here between the PHOBOS 
measurements and LUCIFER simulations. Reasonable descriptions of both the observed DN/Dr) spectra in Au + Au 
as well as the ratios of p± distributions in D + Au and Au + Au at 200 GeV were already obtained with LUCIFER 
simulation. [ToL ITlT] . 

We previously cited [ll[ the theoretical difficulties that arise in producing sufficient flow in a straight-forward 
parton cascade models, e.g. the studies done by Molnar pjj]. This author suggested the possible use of parton 
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coalescence as a mechanism to produce the otherwise too large observed ellipticity. Otherwise, as Molnar notes, the 
parton-parton cross-section must be made as large as 49 mb [l3[ to yield the measured ellipticities. In the sense 
that we are in fact exploiting the properties of a strongly interacting pre-hadronic fluid at a second, but still early, 
stage in our simulation, this hypothesis of Molnar [jjj and our approach do have considerable overlap, practically 
speaking. A relatively large size for the interacting objects present in the cascade dynamics is surely a key ingredient 
in reproducing the required elliptical flow. 

We note from the outset, though, that an appreciable degree of flow, perhaps more than 30%, is already generated 
by the phase I nucleon-nucleon interactions, and subsequently transmitted to the produced pre-hadrons colliding in 
Phase II of the simulation. These initial interactions are mediated by the even larger baryon-baryon cross-sections, 
and the flow in this initial phase is of course an imprint of the geometrical asymmetry which is inherent in ion-ion 
collisions occurring at non-zero impact parameter. 

The approach presented here may be of less validity for the lowest p± charged particles, but becomes increasingly 
accurate with rising p± . Nevertheless we display results from the modeling for even low p± . 

In what follows we outline the dynamics behind our earlier calculations which yielded a successful description of both 
D + Au and Au + Au data at s/s = 200 GeV. Of course the same dynamics is employed in the present determination 
of the ellipticity. Subsequently, actual calculations of elliptical asymmetry are displayed for a variety of kinematic 
cuts in centrality and transverse momentum and compared, in particular, to PHOBOS measurements. 



II. THE SIMULATION 



The simulation code LUCIFER, developed for high energy heavy-ion collisions has previously been applied to both 
SPS energies y/s = (17.2,20) A GeV 0,111 and to RHIC energies ^fs = (56, 130,200)A GeV 0, GH 111 • Although 
nominally intended for dealing with soft, low interaction it is possible to introduce high p± hadron spectra via 
the NN inputs, which form the building blocks of the simulation, and to then examine the effect of re-scattering, and 
concomitant energy loss, on such spectra 0, [111 ]. The simulation is divided into two phases I and II, with most re- 
scattering and explicit two body energy loss restricted to the second and considerably reduced energy stage. The first 
stage sets up the participants, both mesonic and baryonic, their four momenta and positions for the commencement of 
the cascade in II, and hence then also involves energy loss through transfer to the produced particles. The second stage 
energy loss and interactions within a pre-hadronic fluid played a key role in the reproducing the observed suppression 
of the p± distribution in our description of Au + Au collisions. 

The purpose of describing such high energy collisions without introducing the parton structure of hadrons, at least 
for soft processes, is to set a baseline for judging whether deviations from the simulation measured in experiments 
existed and could then signal interesting phenomena. The dividing line between soft and hard processes, the latter 
being in principle described by perturbative QCD, is not necessarily easy to identify in heavy ion data, although many 
authors believe they have accomplished this within a gluon-saturation picture [13, [H, HH • 



A. Pre-hadrons: Production and Hadronisation Time Scales 



Such a mechanism for the suppression of high transverse momentum jets may seem to run counter to the conventional 
pQCD description. There does exist, however, some theoretical justification for a picture in which colourless pre- 
hadrons, may be produced at rather early times in a AA collision, and then play a key role in the further development 
of the system. First, there is the work of Shuryak and Zahed [20( on the persistence of hadron-like states above 
the critical temperature in a dense quark-gluon medium as well as similar results from lattice-driven studies on the 
persistence of the J/'F and other special hadronic states 2l| . 

Recently a much more transparent treatment has been given, on which we now mainly rely, by Kopeliovich, Nemchik 
and Schmidt [22[ as well as Berger [23[ , who directly consider the temporal dynamics of pre-hadron production from 
a pQCD point of view. In particular Kopeliovich et ai, outline the fate of leading hadrons from jets, in deep inelastic 
scattering (DIS) on massive ions and discuss the relevance to relativistic heavy ion collisions. These works have the 
distinct advantage of relying on pQCD. We attempt to reproduce their development as it is germane to our treatment 
of NN, NA and AA interactions. 

Figure (1), essentially borrowed from Reference(22), displays the production of a jet, here at least a moderately 
high energy quark jet, originating in a pp or pA event, as well as the time scales relevant to the process. Kopeliovich 
et al. begin by explicitly considering eP or eA, in which case the initiating particle in the jet creation is an off-shell 
photon (7*). In pp or AA one could equally well employ a gluon at the first vertex. The kinematics [22| remains 
essentially the same then for nuclear induced events and we elaborate this case somewhat. 
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The essential feature evident in the figure is the production of a colour neutral pre-hadron after a time t p , from 
an initially perturbative high p± quark, by what one could label as coalescence of the leading quark with an ambient 
anti-quark. This coalescence incidentally is most probable for co-moving anti-quarks, the scarcity of co-movers with 
increasing p± thus explaining the rapid fall off of even the pp meson spectrum with increasing p± . Also prominent 
in this diagram are the early perturbative gluon radiation. These radiated gluons may in due course create further 
pre-hadrons, and of course there is then concomitant energy loss from the initial quark. Further, for leading partons, 
one expects that the hadronisation time th is considerably longer than t p , involving soft processes that eventually put 
the pre-hadron final on the hadronic mass shell. Octet g<fs might also arise but will not persist or eventually coalesce 
given the repulsive or non-confining forces in such systems. 

One can summarise the kinematic arguments [22| most easily in the rest frame of the struck proton or nucleus. We 
use these arguments to estimate the production time t p for the pre-hadron, which will be a considerably shorter time 
T p in the colliding frame of an AA system. For DIS on a nucleon, on its own or within a nucleus A, the production 
time in Figure (1) is estimated to be 



[E q ] 
[dE q /dz 



(1 - *h), 



(2) 



where Zh is the fraction of quark energy imparted to the hadron. I nteg rating the gluon radiation spectrum one obtains 
for the energy loss per unit length z, a time independent rate [H, H3 rising quadratically with the hard scale Q, 



^ - ^fl>. (3) 



The colour neutralisation time is then given by: 



E, 



*p~qK 1 -^). ( 4 ) 

where Zh = Eh/E q is the fraction of energy imparted to the hadron. 

The hadronisation time is related to the QCD scale factor and is usually estimated as: 

il QCD 

and is much longer than the colour neutralisation time, given that the QCD scale is close to 200 GeV. 

The authors of Reference (22) argue that the brevity of the de-colourisation time scale is a consequence of the 
effects of energy conservation, coherence and Sudakov suppression. 

Importantly, a much reduced time scale for pre-hadron creation is likely to remain true even for non-leading partons 
provided Equation (4) remains valid, i. e. when Q is comparable to E q . Thus many of the radiated hard or moderately 
hard gluons, in the early stages of a pp or AA collision will initiate similar pre-hadrons, and in a nuclear medium 
such large sized objects will be numerous and of critical importance to the dynamics. One should keep in mind that 
the overall multiplicity is not large at p± > 1 GeV/c, where the NN p± spectrum, as seen in Figure (2) [25(, has 
already dropped well below that at the softest transverse momenta measured: by close to two orders of magnitude at 
p± ~ lGeV and for 200 GeV Au+Au as seen in Figure (3) [26| by somewhat more. 

One concludes that the production and hadronisation processes are to some degree separate: with generally t p less 
than tf and frequently much less. Corresponding time scales in, say, the center of momentum frame for an A+A 
system will of course be considerably contracted, given respectively by t p j ~ [r -1 ]t p ,/- Although the colourless 
pre-hadron in Figure (1) is generated by an initially perturbative process with an anti-quark, it is the subsequent 
interactions with other such pre-hadrons that leads to the observed suppression for mesons of appreciable and even 
moderate transverse momentum. The pre-hadron perturbatively begins life with the qq at small relative distance, 
and thus has a small initial size, but in light of the ambient momenta for such high or even moderately high energy 
partons, the transverse diameter rapidly increases to the scale of a typical hadron. Reference (22) suggests the entire 
growth to pre-hadronic size occurs quickly. 

One can pursue the evolution of the system of pre-hadrons via a Glauber-theory [13, H3| based treatment of their 
interactions, or, as we do, via a standard cascade model. The end result is little different: what has been labeled jet 
suppression results. In Glauber theory the hadron-sized cross-sections produce strong absorption: hard pre-hadrons 
simply do not remain in the final state, they are too often absorbed. In the cascade described hereafter the pre- 
hadronic medium is sharply cooled by the interactions and instead of one hard meson, many softer mesons appear at 
lower p± . 
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To obtain the final cascade yields for both D+Au Q and Au+Au [Tl|, it's essential that the characteristic time t p 
is considerably less than tf. Indeed, this constraint, and the appropriately large pre-hadron interaction cross-sections 
also are critical in generating the surprisingly large elliptical flow that has been measured at RHIC [J, 0, H 
Kopeliovich et al. provide a perturbative QCD justification for both conditions: early de-colourisation and large 
pre-hadron interaction cross-sections. 



B. Simulation Dynamics 



For completeness we present a brief overview of the dynamics of our Monte Carlo simulation, which has in fact 
already been described extensively in earlier works 0, [ill , fill ] . Many other simulations of heavy ion collisions exist 
and these are frequently hybrid in nature, using say string models in the initial state [H, \2^L 130. l3ll 32l. 33l. 34l. f35l| 
together with final state hadronic collisions, while some are either purely or partly partonic 
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nature. Our approach is closest in spirit to that of RQMD and K. Gallmeister, C. Greiner, and Z. Xu [42| as well as 
work by W. Cassing [43| and G. Wolschin [44J. The latter authors seek to separate initial perhaps parton dominated 
processes from hadronic interactions occurring at some intermediate but not necessarily late time, and/or rely on 
some hadronic dynamics. 

Whereas in the simulations of the D + Au and Au + Au ion-ion collisions our purpose was to present a background 
hadronic setting from which deviations in observations would signal interesting QCD phenomena, the measured flow, 
which is certainly appreciable, seems ab initio to demand a pre-hadronic source. In stage I in the incoming target 
and projectile nucleons interactions are tracked, while in II produced particles, considered as pre-hadrons interact 
in a hadron-like cascade. Basic inputs in the entire simulation [14L I la . llQ] are the measured hadron-hadron cross- 
sections [H, |H, HgJ as functions of energy and the production multiplicities taken to conform to KNO [13] scaling. 
The elementary NN dynamics is split, as in the classic work of Goulianos on the phenomenological treatment of 
NN (48|, into elastic, single diffractive (SD) and non-single diffractive NSD). 

No energy loss is permitted in I, except at the completion of this stage when production of pre-hadrons is initiated, 
the collisions occurring essentially parallel in time. A history of each collision is retained, from which one can construct 
the particle p± from random walk and in the present case can also infer, for non-vanishing impact parameter, an 
asymmetry in the transverse momentum distributions. This p± and attached asymmetry can be appropriately passed 
on to the pre-hadrons generated at the end of phase I. 

The transition from stage I to stage II is accomplished with the guidance of the two-body inputs. The multiplic- 
ities, energies and character of the produced particles are determined, as well as the initial, transverse-asymmetric 
momentum distribution. The entities produced in I are principally a set of vector pre-hadrons, and also after a time 
Tf decay into the normal "stable" mesons. The decays take account of the vector nature of the pre-hadrons. In this 
work the baryons are limited to the usual flavour octet, the comparison with elementary NN measurements severely 
limits any freedom in the mass, number and spectra of the produced prcmesons. 

Pre-hadrons, which when mesonic presumably consist of a spatially loosely correlated quark and anti-quark pair, 
are given a mass spectrum between m ~ m^ and m ~ 1.2 GeV, with correspondingly higher upper and lower limits 
allowed for pre-hadrons including strange quarks. The Monte-Carlo selection of masses is then governed by a Gaussian 
distribution, 



P(m) = exp(— (m — mo) 2 /w 2 ), (6) 

with mo a selected center for the pre-hadron mass distribution and w = mo/4 the width. For non-strange mesonic 
pre-hadrons mo ~ 700 MeV, and for strange mo ~ 950 MeV. Small changes in mo and w have little effect since the 
code is constrained to fit hadron-hadron data. Too high an upper limit for mo would destroy the soft nature expected 
for most pre-hadron interactions when they finally decay into 'stable' mesons. 

To reiterate: the cross-sections in pre-hadronic collisions are taken to be the same magnitude as hadronic cross- 
sections, at the same center of mass energy, thus introducing no additional free parameters into the model. Where 
hadronic cross-sections or their energy dependences are inadequately known we employ straightforward quark counting 
to determine the cross-section magnitudes. Further details on this early stage of the simulation can be found in earlier 
work & El Eft 

It should also be noted that the rather small DN/Dr] observed by PHOBOS [49| was in fact predicted in 
our simulation directly from the rise in multiplicities with energy inserted into LUCIFER via the experimentally 
constrained NN inputs. This very small rise puts some doubts on the achievement of a standard QCD plasma at 
RHIC. The increase of entropy expected in a phase transition from confined to unconfined partons should show up as 
a rather sharp increase in produced mesons. 
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III. STAGE II: FINAL STATE CASCADE 

Stage II is a straightforward two body cascade in which the pre-hadrons, and of course any normal hadrons present, 
may interact and decay. The pre-hadron decay time, taken at rest or in fact in the colliding Au+Au frame uniformly 
as Tf ~ 1.0 fm., can be viewed as a hadronisation or formation time. Appreciable energy having being finally 
transferred to the produced particles these 'final state' interactions occur at considerably lower energy than the initial 
nucleon-nucleon collisions in stage I. The final cascade of course conserves energy-momentum and leads to additional 
transverse energy and meson production. For Au+Au, the effect of pre-hadron interactions is appreciable, greatly 
increasing multiplicities and total transverse energy, E± , through both production and eventual decay into the stable 
meson species. These 'final state' interactions of stage II were the principal agent in the suppression of high p± 
production [{j ITTj|. 

The spatial positions of the particles at this time could be assigned in various ways. We have chosen to place 
the pre-hadrons inside a cylinder, extending in the transverse and longitudinal directions far enough to contain the 
final stage I positions of the initial baryons in the group. We then allow the cylinder to evolve freely longitudinally 
according to the corresponding momentum distributions, for a fixed time r p , i. e. roughly the pre-hadron production 
time in the colliding frame. At the end of r p the total multiplicity of pre-hadrons is limited so that, given normal 
hadronic sizes appropriate to meson-meson cross-sections ~ (2/3)(47r/3)(0.6) 3 fm 3 , hadrons do not overlap within the 
cylinder. Such a limitation in density is consonant with the notion that produced hadrons can only exist as particles 
when separated from the interaction region in which they are generated [50j. One may conclude from this that the 
pre-hadronic matter acts like, at its creation, as an incompressible fluid, viz. a liquid, as depicted in the earliest 
calculations with LUCIFER [1, G3 . 

IV. FLOW DETERMINATION 

We have run simulations for central, mid-central and peripheral cuts used by PHOBOS Q to correspond to cen- 
tralizes of 3%-15%, 15%-25%, and 25%-50% respectively. Calculated results are displayed in Figure(4) alongside the 
PHOBOS determinations. The theoretical curves shown are best smooth fits to the LUCIFER results, which of course 
display statistical errors. Given the larger non-statistical experimental errors and the difficulty of calculation at large 
one probably cannot yet take very seriously the apparent upturn in v-i at the extreme of the pseudo-rapidity range. 
The agreement between measurement and simulation is quite striking. With the code already completely specified 
and constrained by the previous calculations of rapidity and transverse momenta spectra [H, [Til ] for both D+Au and 
Au+Au, there is nothing to adjust to yield the asymmetric elliptical flow. 

Of course we must also confront the dependence of ellipticity on transverse momentum, pj_, evidenced in the 
experimental data. This comparison is given in Figure(5). The present model, without intervention, handles this 
behaviour more than adequately. There is an agreement in both magnitude and p±_ variation, in particular the clear 
saturation of ellipticity above 1-2 GeV/c. For the moment we cannot follow the functional variation of V2 to very high 
p±, since it requires a considerable increase in computation time. Nevertheless, the deviation of our model from a 
naive hydrodynamic picture is already evident. More importantly there is unambiguous confirmation of the presence 
of hadron-like degrees of freedom, at an early stage in the 200 GeV Au+Au collision. The observed magnitude of flow 
could not have obtained without hadron sized interaction cross-sections in both phase I and II of the simulation. 

The non-hydrodynamic saturation of V2 with increasing p^ presents an interesting puzzle. Our best guess is 
that at the higher transverse momenta there is less time for interaction within the stage II hadronic medium and 
correspondingly less tying of the ellipticity to the initial geometry for the highest p± mesons. Also decays of pre- 
hadrons play an important role here: both decreasing the flow at low p± and increasing it at high p± . 

There are two clear lessons to be learned. First that the second phase, the pre-hadronic cascade, must begin early, 
at a time t p ~ (T)~ 1 t p with T the average Lorentz factor for pre-hadrons in the collision center of velocity frame for 
the colliding gold nuclei. This is in accord with what was found earlier in simulations of the single particle spectra in 
d+Au and Au+Au. Second that large cross-sections must be employed in phase II and indeed also in the initiating 
NN collisions in phase I. 

The flow is apparently a consequence of and requires large interaction cross-sections, between objects which cascade 
starting at early times. This picture either obviates or is in accord with the phrase "strongly interacting quark-gluon 
plasma" which has been popularly employed to described the dense matter produced in the RHIC energy ion-ion 
collisions, depending on one's point of view. The aptness of the description of the pre-hadron medium as a liquid was 
made clear in our earlier work, where the pre-hadron multiplicity was limited by a constraint that the pre-hadrons 
could not overlap d, [HI EH • 
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A. Time and Interaction Dependence of Ellipticity. 

By way of a theoretical experiment, we have in Figure (6) explored the variation in produced flow with the initial 
delay r p before the start of phase II and also with the magnitude of the meson-meson cross-section in this phase II. The 
results are clear. Halving the cross-sections decidedly reduces the flow, v 2 , to a value well below the measured results. 
A similar effect occurs when the pre-hadron production time t p is increased. In a dense many body system there 
is a natural screening with the cross-section being effectively limited by that corresponding to the average distance 
between particles. Should the production of flow occur early on, when the density is high, the limiting cross-section 
will be ~ ira 2 with a the average interparticle distance. As a result little effect will be evident from a cross-section 
increase. 

The t p dependence exhibited in Figure(6) is striking, the flow diminishes very appreciably with a large delay in 
the commencement of phase II. The geometrical, free expansion of the system reduces both the densities and the 
interaction frequencies, thus increasing the mean free path for pre-hadrons; in hydrodynamic terms this lessens the 
pressure which produces the flow. The connection between ellipticity and the initial spatial configuration of the 
ion-ion collision is thus lost by the time the interactions begin. One should note that the cascade calculations in an 
appropriate limit should yield hydrodynamics. Of course the hydronamic limit for flow is achieved experimentally only 
for low pj_ and in Figure(5) both measurements and LUCIFER simulation drop well below the earliest hydrodynamic 
predictions at higher low p± ~ 2 GeV/c. 

Very soft early gluons and quarks are not a part of the colour neutral medium, but their eventual hadronisation is 
approximated by our treatment of soft hadronic processes, guided by known hadron-hadron interactions and produc- 
tion multiplicities. Finally, completely eliminating phase II reduces the ellipticity v 2 by some 25-30%, depending on 
the centrality, the remaining flow being ascribable to the initial nucleon-nucleon interaction which already occurred 
in phase I. 

V. COMMENTS 

Overall the very natural fashion in which the elliptical flow is produced in the model speaks strongly to the 
important role played by pre-hadrons, i. e. colour neutral objects having hadron-sized interaction cross-sections, and 
their presence from an early time in Au+Au collisions at RHIC energies. Elliptical flow is, after all, the one truly 
collective phenomenon to have been observed at RHIC, and it is precisely such collective properties that may indicate 
the existence of some medium or fluid, in this case a plasma apparently characterized by objects having hadron sized 
interaction cross-sections during much of its existence. The very earliest times in an A+A collision, tracked in the 
first stage of our simulation, also play a non-negligible role in the production of elliptical flow. The nucleon-nucleon 
interactions occurring in this stage likely lead to a Cronin-like [5l| enhancement [ll| of the high p± spectrum and 
also transmit the spatial geometry of the collision to the produced pre-hadrons. It would be interesting to evolve this 
first stage with allowances made for mutual crossing between parallel quark-gluon and hadronic simulations. Future 
work will pursue such an avenue. Also, in agremeent with the PHOBOS observations Q, the simulation produces 
only very small direct flow, i. e. v\, at s 1 ' 2 = 200 GeV. 
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FIG. 1: Schematic drawing, borrowed from Kopeliovich et al. [23], of the perturbative formation of a pre-hadron from an 
off-shell gluon g* incident on a quark in the rest frame of nucleus B in an A+B or P+B event. The quark radiates gluons and 
eventually combines with a perturbatively produced anti-quark to form at first a small colourless pre-hadron which rapidly 
expands to hadronic size. The time t p signifies the production time of the pre-hadron and th its later hadronisation time. 
It is argued in the text that in general t p is much less than th, and the existence of two such time scales is critical to the 
parton-hadron dynamics. 




FIG. 2: pp Pseudo-rapidity spectra: Comparison of UA1 minimum bias 200 GeV NSD data [25( with an appropriate LUCIFER 
simulation. The latter is properly constrained by experiment and is an input to the ensuing AA collisions; thus does not 
constitute a 'set' of free parameters. 




FIG. 3: Central PHENIX ir 200 GeV for Au+Au vs simulation. Curves for different choices of the production time t p differ 
very little, since in effect the cascade effectively begins somewhat later, near 0.25 — 0.35 fm/c and continues much longer to 
tens of fm/c. Centrality for PHENIX is here 0% — 10%, roughly for impact parameters 6 < 4.25 fm. in the simulation. 
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FIG. 4: PHOBOS 200 GeV charged flow, V2, as a function of pseudo-rapidity for peripheral, mid-central, and central versus 
LUCIFER simulations, with actual points calculated indicated in central and peripheral but for an already busy graph only a 
fit is shown for mid-central. 
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FIG. 5: PHOBOS 200 GeV V2 as a function of transverse momentum for 0-50% centrality versus LUCIFER simulations for 
both charged mesons and baryons. Hydro calculations are also indicated 
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FIG. 6: Variation in LUCIFER central simulations with scale changes for the overall pre-hadron-pre-hadron cross section by a 
factor of one-half and the time factor by 2.0. This decrease in interaction strength, especially, and lengthening of the time lead 
to considerable reduction in the elliptical flow. 



